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Abstract. We report the results of a new large scale sur- 
vey of the Orion-Monoceros complex of molecular clouds 
made in the J = 1 — > line of 12 CO with the Harvard- 
Smithsonian 1.2 m millimetre-wave telescope. The survey 
consists of 52 288 uniformly spaced spectra that cover an 
area of 432 deg 2 on the sky and represent the most sen- 
sitive large-scale survey of the region to date. Distances 
to the constituent molecular clouds of the complex, es- 
timated from an analysis of foreground and background 
stars, have provided information on the three dimensional 
structure of the entire complex. 
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1. Introduction 

Orion-Monoceros is one of the most important and best 
studied regions of star-formation. Much of its importance 
derives from its favourable position on the sky and its 
proximity: it is located ~ 15° below the Galactic plane 
and toward the outer Galaxy — well away from the poten- 
tially confusing molecular gas of the plane and the Galac- 
tic centre; at a distance of ~ 450 pc, it is relatively close, 
and its constituent clouds can be observed with good lin- 
ear resolution even with a radio telescope of modest size. 
The radial velocities, proper motions and parallaxes of the 
stars associated with the molecular gas are observable, 
and therefore the distances to the clouds and the three- 
dimensional structure of the complex can be estimated. 

Orion-Monoceros provides an excellent laboratory for 
studying the interaction between massive stars and the 
interstellar medium (ISM). Over the course of the last 
12Myr the molecular clouds in the region have been 
shaped, compressed and disrupted by the powerful ion- 
izing radiation, stellar winds and supernova (SN) explo- 
sions of the young massive stars in the Orion OB associ- 
ation (Cowie, Songailia & York 1979; Bally et al. 1987). 
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Compression of the clouds may have triggered bursts of 
star-formation (Elmegreen & Lada 1977) that are traced 
by the four semi-distinct Ori OBI subgroups. The old- 
est subgroup, OB la, lies in a region about 10° from the 
main clouds long-ago cleared of molecular gas, while the 
youngest OB Id is currently forming the stars in the Orion 
nebulae on the near side of the Orion A cloud. As well 
as disrupting the molecular gas in Orion-Monoceros, the 
stars of Ori OBI have also inflated a huge elongated bub- 
ble of hot X-ray emitting gas, the Orion-Eridanus super- 
shell (Heiles, Haffner & Reynolds 1999), which is centered 
on the Orion star-forming region 15°-20° (~ 150 pc) be- 
low the Galactic Plane. As the bubble has expanded it has 
swept up a massive shell of dust and gas that has become 
an optically bright crescent (Barnard's Loop) as it moved 
into the steep density and pressure gradients associated 
with the Galactic gas layer. Fig.^shows a finder diagram 
for Orion-Monoceros. The colourscale represents the in- 
tegrated Ha emission as measured by the Wisconsin Ha 
mapper (WHAM - Reynolds et al. 1998), and the con- 
tours show the boundaries of the main molecular clouds 
as traced by the CO survey presented here. The positions 
of the O and B stars from the four Ori-OB 1 subgroups, 
from Brown et al. (1995), are also plotted. 

The first observations of molecular emission from 
Orion-Monoceros, which were also the first detections of 
interstellar CO, were made toward the Orion nebula by 
Wilson, Jefferts & Penzias (1970) using the J = 1 — > 
line of the normal isotopic species. Subsequent studies us- 
ing this line by Kutner et al. (1977), Chin (1977), and 
in particular Maddalena et al. (1986) were conducted pri- 
marily to determine the spatial extent of the complex. It 
was found that the complex contains three giant (10 5 M Q ) 
molecular clouds (GMCs), Orion A, Orion B and Mon R2, 
two long filaments which extend approximately 10° from 
the cloud complex to the Galactic plane, and numerous 
smaller clouds. Parts of the complex have subsequently 
been observed in CO at higher resolution either using the 
J = 1 — ► line or other isotopes or transitions. Most 
importantly, the majority of Mon R2 was mapped in CO 
(J = 1 -> 0) by Xie & Goldsmith (1994), the central 
part of Orion A was mapped using 13 CO ( J = 1 — > 0) by 
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Fig. 1. Finder diagram for Orion-Monoceros. The colourscale represents the integrated intensity of ionized gas as 
traced by the WHAM Ha survey (Reynolds et al. 1998). The bright crescent of emission is Barnard's loop which 
delineates the optical boundary of the Orion-Eridanus superbubble. The positions of O and B stars from the Ori OBI 
association (Brown et al. 1995) are plotted with symbols the sizes of which are proportional to the luminosity of the 
stars. The four different symbols correspond to the four Ori OB 1 subgroups. The boundary of the Orion-Monoceros 
molecular complex is marked by the 2.06Kkms~ 1 (3a) contour of the new CO survey. 



Bally et al. (1987), by Nagahama et al. (1998), and in the 
J=1^0andJ = 2^1 lines of 12 CO, 13 CO and C 18 
by Castets et al. (1990). In addition, the southern part of 
Orion B was mapped in CO and 13 CO ( J = 2 — » 1 and 
J = 3 — > 2) by Kramer, Stutzki & Winnewisser (1996), 
and the Orion A and B clouds were mapped using the 
J = 2 — > 1 line of CO (Sakamoto et al. 1994). However, to 
date the only survey which has covered the whole region 
has been that of Maddalena et al. (1986), a survey which 
is non-uniform in resolution, sensitivity and sampling, and 
is heavily undersampled in many places. 

This paper presents a new CO (J = 1 — > 0) survey 
of the Orion-Monoceros complex which consists of 52 288 
spectra uniformly spaced in Galactic longitude and lat- 
itude. It is generally four times more sensitive than the 



Maddalena et al. (1986) survey, since it contains four times 
as many spectra, each with lower rms noise. The new ob- 
servations cover a total area of more than 400sqdeg on 
the sky and map all of the major clouds of the Orion- 
Monoceros complex with the exception of the A-Orionis 
ring, which was the subject of a separate study with the 
1.2 m telescope (Lang et al. 2000). 

This paper is organised as follows: In f|3]the observa- 
tions and data reduction procedure are described. In 
the methods used to estimate the masses of and the dis- 
tances to specific regions in the Orion-Monoceros complex 
are outlined. In ^1 these specific regions are described in 
detail. §5\ presents an overview of the large scale struc- 
ture of the complex and describes a possible formation 
scenario. A summary and conclusions are presented in ^1 



B. A. Wilson et al.: A Uniform CO Survey of the Molecular Clouds in Orion and Monoceros. 



3 




I I I i i I I I I I I I i i I I I r I L. I I i i I - 

225" 2211* 21&= 210" 205" 200 

Galectic Longitude 



Fig. 2. Log velocity-integrated ( J = 1 — > 0) CO map of Orion-Monoceros. The integration range is —3.0 < Ulsr < 
19.5kms _1 . Contours start at 0.313 (3cr) and go up in steps of 0.313Log(Kkms _1 ). The dashed line marks the limit of 
the full-resolution survey. The region on the high longitude side of the dashed line has been sampled with beamwidth 
spaced observations; that on the low longitude side has been observed twice with the second set of observations offset 
by half a beamwidth in longitude and latitude. 



2. Observations 

The 1.2 m millimetre wave telescope at the Harvard- 
Smithsonian Center for Astrophysics was used for all the 
J = 1 — > 12 CO observations presented here. This tele- 
scope has a beamwidth (FWHM) of 8'4 and a velocity 
resolution of 0.65 km -1 at 115.2712GHz, provided by a 
256 channel filter bank. Spectral line intensities were cali- 
brated against an ambient temperature blackbody, a stan- 
dard Eccosorb chopper wheel (Davis et al 1973), which was 
rotated in front of the feed horn at 20 Hz for 1 s before each 
observation. All line intensities reported here are in units 
of main beam brightness temperature T m b = T^/77, where 
T^ is the chopper calibrated antenna temperature (Pen- 
zias et al. 1973) (i.e. the antenna temperature corrected 
for atmospheric attenuation, ohmic losses and rearward 



spillover), and rj is the main beam efficiency, 0.82. The 
emission in Orion is so strong and extensive that Tmb 
will at some locations overestimate the intensity in the 
main beam by a factor up to 5% owing to emission in the 
near sidelobes. The detector was a liquid- helium cooled 
SIS receiver with a single sideband noise temperature of 
- 75 K (Pan 1983). 

Measurements of the sky brightness as a function of 
elevation were obtained by antenna tipping at least every 
six hours and more frequently during periods of change- 
able weather, and were fitted to the two layer atmospheric 
model of Kutner (1978) to determine the temperature and 
opacity of atmospheric water vapour. Further details on 
the receivers, observational technique and calibration of 
the 1.2 m telescope can be found in Dame et al. (1993). 
During the course of the observations total system tern- 
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Fig. 3. Channels maps. Each map shows the integrated emission from five velocity channels (total width 3.25 kms ) 
from —1.0 to 18.5kms _1 . Log contour levels start at 0.1 (3er) and go up in steps of 0.3. 



peratures referred to above the atmosphere ranged from ~ Frequency switched observations were carried out over 

900 K at the lowest elevations observed (30°) to ~ 500 K four consecutive observing seasons from February 1999 

at the highest (~ 60°). Integration times were automat- to April 2002, using two 0?125 grids that were offset by 

ically adjusted to achieve an rms noise temperature per 0?0625 in Galactic Longitude and Latitude. The beam 

channel of ~ 0.26 K. spacing of the final composite grid is 0?09, and the x and 
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Fig. 4. Emission weighted mean velocity map of Orion- 
Monoceros. This figure clearly shows the velocity gradi- 
ents along Orion A and the Southern Filament and high- 
lights regions with unusual kinematics such as the Scissors 
and NGC 2149. 



y axes of the grid are tilted at 45° with respect to I and b. 
The frequency switched spectra were folded and 5 th order 
polynomial baselines fitted to the portion of each spectrum 
that was found to be free of emission. For each observing 
session of ~ 6 hours a model of the CO emission feature 
associated with the mesospheric layer of the Earth's atmo- 
sphere was fitted to all spectra that were free of Galactic 
CO emission (Lang 1997). This model was then used to 
remove the mesospheric CO line from all the spectra taken 
for that session. Finally the oversampled observations were 
smoothed by convolution with a two-dimensional Gaus- 
sian the size of the telescope's main beam in order to lower 
the final rms noise to ~ 0.14 K. The parameters of the pre- 
smoothed observations are summarized in Tabled 

The CO observations are summarized in Fig. [3 which 
shows the spatial distribution of Wco (Kkms -1 ), the 
velocity-integrated line intensity. The i>lsr integration in- 
terval of —3 to 19.5 km s" 1 was chosen after an examina- 
tion of the channel maps in order to include all the regions 
in the spectra where signal was present. The higher sen- 
sitivity per unit area of the new survey is apparent in the 
amount of emission that is now detected from the areas 
between the main molecular clouds. Although weak, this 
emission bridges the clouds spatially and without disconti- 
nuity in velocity, indicating that none of the major clouds 
are entirely isolated. The velocity structure of the region 
is shown in the channel maps, Fig. |21 and the emission- 
weighted mean velocity map Fig. 0] These two figures in- 
dicate that the majority of the emission is confined to the 
relatively narrow velocity range 2 — > 15kms _1 ; although 
individual clouds such as Orion A and the Southern Fila- 
ment have systematic velocity gradients, there is none for 



Number of spectra 
Galactic longitude 
Galactic latitude 
Angular resolution 
Sampling interval 
Velocity resolution 
Sensitivity 

Frequency switching interval 



52 288 
200°^226° 
-22° ^-3° 
8.'4 
0?088 

0.65 kms -1 
T mb = 0.26 K rms 
15 MHz 



the complex as a whole. The structure and kinematics of 
the individual regions within the complex are discussed in 
detail in 

3. Molecular Masses and Cloud Distances 

In spite of clear evidence of saturation in the denser cores 
of molecular clouds, the 12 CO line remains the best over- 
all tracer of mass in giant molecular clouds, since most 
of the mass is contained in the lower-density envelopes 
(Cambresy 1999). Cloud masses in the Orion-Monoceros 
complex were estimated using the standard assumption of 
a linear relationship between the velocity integrated CO 
intensity, Wco , and the molecular hydrogen column den- 
sity, iV(H 2 ), taking 

N(R 2 )/W C o = (1-8 ± 0.3) x 10 20 cm" 2 KT 1 km" 1 s" 1 

as derived by Dame et al.(2001). This yields 

Mco/M© = 1200S CO < C 

where d is the distance to the cloud in kpc, and Sco is the 
CO emission integrated over velocity and the angular ex- 
tent of the cloud in Kkms" 1 arcdegrees 2 . A mean atomic 
weight of 1.36 is assumed in order to account for helium 
and heavier elements (Allen 1973). 

Historically the distances to Orion A and B have been 
assumed to be the same as the distance to the Orion neb- 
ula which has been estimated photometrically as 435 pc 
(Warren et al. 1977), and as 480±80pc via a study of 
the radial velocities and proper motions of masers in the 
OMC1 core (Genzel et al. 1981) 1 . Although these esti- 
mates may be accurate, the implicit assumption that all 
the clouds are at the same distance is unsatisfactory. At 
450 pc the Orion A cloud would be over 100 pc long, and 
although the region near to the Orion nebula is likely to be 
at a distance close to the photometric estimate, the other 
end of the cloud may be significantly nearer or further 
away. 

In order to address this problem the parallaxes of stars 
observed with the astrometric satellite Hipparcos have 
been used to estimate distances to the molecular clouds in 



1 Mon R2 is generally believed to be further away at 
830±50pc (Racine 1968; Herbst & Racine 1976) and not di- 
rectly related to the other clouds in the complex. 
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Table 2. Masses and distances for the molecular clouds in the Orion-Monoceros complex. 



Region 


Distance 


Previous 
Estimate 


Reference 


Molecular 
Mass x 10 3 




(pc) 


(PC) 




(M ) 


Orion A: 








105.1 


1 


521+if 






12.3 


2 


465j£ 


480 ±80 


Genzel et al. (1981) b 


69.5 


3 


412±^ 4 






23.3 


NGC 2149 


425 c 






12.2 


Orion B: 








82.3 


1 near 


422±g 


400 


Anthony -Twarog (1982) 


11.9 


1 far 


514±i° 6 






53.4 


2 


387±7i 2 






18.1 


Scissors 


- 150 rf 






0.340 


Orion East 


l2Q±f 






0.129 


Monoceros R2 


80011™ 


830 ±50 


Herbst & Racine (1976) e 


106.2 


Southern Filament: 








17.4 


1 (Crossbones) . . . 


465+T 








2 


457+if 








Northern Filament . . . 


393+^8 






17.5 



a la uncertainties 

b Also ~ 435 pc Warren & Hesser (1977) 

c No background stars. Mean of distances to Orion A (3) and Southern Filament (1) 
adopted. 



No foreground stars so distance estimated from background stars only. 
e Also - 950 pc Racine k van den Bergh (1970) 



Orion-Monoceros. Individual stars are identified as being 
in front or behind a particular molecular cloud by com- 
paring its colour with two theoretical models: the colour 
of the star if it is reddened by dust associated with atomic 
and molecular gas (i.e it is in the background of the cloud), 
and the reddening associated with just atomic gas (i.e. it is 
in the foreground of the cloud) . The parallaxes of the stars 
that can be assigned unambiguously to the foreground or 
the background of a molecular cloud are then used to con- 
strain the distance to that cloud. 

In general Hipparcos parallaxes for single stars are un- 
reliable for distances greater than ~ 200 pc, but by using 
groups of foreground and background stars and carefully 
modelling the probability distribution function (pdf) for 
the cloud distance it is possible to make robust distance 
estimates out to approximately 500 pc. Furthermore, in 
some cases it has been possible to divide the larger clouds 
into multiple sections and make separate distance esti- 
mates for each part. In these cases the velocity structure 



of the complex was used as a guide when sub-dividing a 
cloud. 

Table [3 lists the distances and masses of the various 
clouds in Orion-Monoceros estimated from Hipparcos par- 
allaxes. Most of these distances agree, within the errors, 
with previous estimates made using other techniques. The 
distance determination technique is described in more de- 
tail together with specific examples for Orion and other 
local molecular clouds in Wilson et al. (2005). 

4. Specific Regions 

4-.1. Orion A 

The Orion A cloud, the largest in the Orion-Monoceros 
complex, is located ^19° below the Galactic plane at a 
distance of ~ 450 pc from the Sun. It has a surface area of 
~ 31.5 deg 2 (~ 2200 pc 2 ) and if cylindrical in shape, has a 
depth of ~ 20 pc. Fig. [3(a) shows the velocity-integrated 
CO emission from Orion A which is somewhat cometary 
in appearance, particularly above the 24 K km s _1 contour. 
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Fig. 5. (a) Wco map of Orion A. The integration range is 0.0 < «lsr < 15.0kms _1 . Contours start at 1.68Kkms _1 
(3a) and subsequent levels are at 8, 13, 23, 33, 43, 63, 83 and 103(7. The dashed horizontal line marks the nominal 
boundary between Orion A and Orion B, the dotted vertical lines separate the three regions for which distances 
were estimated and the cross marks the position of the Orion Nebula, (b) Longitude-velocity map of Orion A. The 
integration range is —22° < b < —17° and the contours start at 0.468 Karcdegree (3a) with subsequent levels at 3, 8, 
13, 18, 28, 38, 48, 63, 78 and 93ct. 



The apparent head of the comet is the compact knot of 
emission at I — 209° and b = — 19?25, and the ridge of 
dense gas that points away from this peak (and also away 
from the centroid of the OB association) , becomes a wider 
lower density tail toward higher longitudes. 

The longitude-velocity map in Fig. (b) shows the 
large-scale velocity gradient that runs along the length 
of the cloud. First reported by Kutner et al. (1977), this 
gradient has been confirmed by all subsequent large-scale 
observations of the cloud. Its origin has been attributed to 
both rotation (Kutner et al. 1977; Maddalena et al. 1986) 
and to large scale expansion driven by the stellar winds 
of the Ori OBI association (Bally 1987). In the first case 
the rotation is about an axis perpendicular to the Galactic 
plane and opposed to the rotation of the Galaxy. Heyer 
et al. (1992) argue that this model is unlikely to be cor- 



rect as the velocity-gradient (as shown by their data) is 
not sufficiently continuous over the length of the major 
axis to be consistent with the cloud's rotational motion. 
Specifically, there exist sharp velocity shifts of ~ 1 kms -1 
between adjacent spectra at many locations. Bally et al 
(1987) and Heyer et al. (1992) also identify many regions 
within Orion A that exhibit very sharp emission gradi- 
ents and velocity shifts which they attribute to expand- 
ing shells driven by HII regions excited by young massive 
stars. However, the velocity gradient as shown in Fig. [3(b) 
extends over a distance of ~ 100 pc and it is unlikely that 
such a coherent structure is the result of purely internal 
processes. It is possible, however, that the gradient was 
formed by the passage through the cloud of a much larger 
shell, similar to those driven by the Ori OB association. 
It is noted that it is not possible to see the huge Orion A 
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outflow in the longitude-velocity (Fig. 0(b)) as the data 
has been integrated over 5°. However, the Orion A outflow 
can be seen clearly in narrower position- velocity slices (see 

Distances to Orion A were determined for each of 
three sections; delineated by the vertical dotted lines on 
Fig. 03(a) at 208° and 213°, using the Hipparcos constrain- 
ing star method described in §3\ The distances were esti- 
mated to be 520 pc, 465 pc and 410 pc for Regions 1, 2 and 
3 respectively. Although adjacent distance estimates are 
roughly equal (c.f. Table EJ, they imply that the cloud is 
oriented in a direction consistent with the rotation model 
of Kutner (1977) and Maddelena (1986). However, this 
model is unable to explain the velocity jump between Re- 
gions 1 and 2. A more likely explanation for the differ- 
ent distances is suggested by the locations of the stars 
of the OB lb subgroup (c.f. Fig. ^) which are believed to 
be located at a distance of ~ 360 pc (Brown et al. 1994) 
in front of the low-longitude end of the cloud. A model 
where Region 1 is accelerated away from the Sun by the 
energy released from the massive stars of OB lb is en- 
tirely consistent with the kinematics of the cloud. Fur- 
thermore there are several small clouds (for example near 
I = 202?5,6 = -17°) located toward the OB lb associ- 
ation that also have high velocities consistent with this 
model. 

If the low-longitude end of Orion A has been acceler- 
ated away from the Sun by the stellar winds of the OB lb 
subgroup then we would expect to see a corresponding 
lateral acceleration of the gas away from the centre of the 
subgroup. The —1.0 — > 8.8 km s -1 channel maps (Fig. 01 
show that the area around I « 205°, b « —19° (the part of 
the cloud closest to the OB lb subgroup) has been cleared 
of molecular gas at low velocities (< 5kms _1 ). In this ve- 
locity range the gas in Orion A, and on the low-latitude 
side of Orion B, may also have been compressed laterally, 
a process that has given rise to the pronounced step in 
the Orion A longitude-velocity map at ~208?5. The ve- 
locity step is coincident with both the peak CO emission 
and the position of what is presently the most active site 
of star-formation within the complex, the Orion nebula 
(Hillenbrand 1997). The fact that this region of high den- 
sity gas (at ~ 465 pc) and the young stars of the Orion 
Nebula (at ~ 435 pc; Warren & Hesser 1977) appear to 
be at similar distances is not entirely unexpected. It is 
difficult to avoid the conclusion that the OB lb subgroup 
compressed the gas in the proto-Orion A cloud which trig- 
gered the formation of the OB lc and OB Id associations. 

Moving along Orion A away from the centre of the 
OBlc subgroup and toward higher longitudes, the gas be- 
comes more diffuse and there is a marked decrease in star- 
formation. Region 3 of the Orion A cloud, the area be- 
tween I w 213° and I « 217°, contains two components 
that are distinct both in position and velocity. Maddalena 
et al. (1986) suggested that the second, higher velocity 
component belonged to a separate group of clouds associ- 



ated with NGC 2149 that is not directly related to Orion 
A. It will be shown in section ED that the NGC 2149 
clouds in fact form a coherent structure with the high lon- 
gitude end of the Orion A cloud and part of the Southern 
Filament. 

4-2. Orion B 

Orion B is the second largest cloud in the Ori-Mon com- 
plex, at least in apparent size. The cloud subtends an area 
of <~ 25.7 deg 2 , which at 450 pc is 1500 pc 2 , and it has a 
mass of 0.8 x 10 5 M . The CO emission from Orion B is 
summarized in Fig. El which shows (a) the Wco map and 
(b) the latitude- velocity map. The dashed line at b = —17° 
on the Wco map marks the nominal boundary between 
the Orion B and Orion A clouds. Emission from the "Scis- 
sors" (discussed below) was not included in the latitude- 
velocity map. 

The emission from Orion B peaks at I = 206? 5, b = 
-16?5, toward the HII region NGC 2024 and the reflec- 
tion nebula NGC 2023. A ridge of high density gas ex- 
tends from the CO peak along the very well defined low- 
latitude edge of the cloud. The low-longitude side of the 
cloud, the side closest to the centres of the OB la and 
OB lb subgroups, also shows a sharp edge in CO emission. 
These sharp edges suggest that the boundary between the 
molecular and the atomic gas is seen close to edge on. By 
contrast there is a network of increasingly diffuse ridges 
and filaments that trail off toward the opposite side of the 
cloud, in directions that point away from the centre of the 
OB association. 

It can be seen from the velocity-latitude plot in Fig. E] 
(b) that the gas in Orion B does not possess a systematic 
velocity gradient and that the cloud as a whole is cen- 
tered at v ~ lOkms . However, the majority of spectra 
from the region below b — —14° contain additional emis- 
sion at a lower velocity (v ss 7kms~ 1 ). This excess can 
also be seen in the low- velocity channel maps (Fig. |3J). 
Many of the CO lines from this region arc cither double 
peaked or unusually wide and lopsided with the excess at 
low velocities. Toward higher Galactic latitudes, around 
I = 204? 5, b = —11? 5, a second cloud is superimposed 
upon Orion B. Associated with the dark nebulae L 1621 
and L 1622 (Lynds 1962) and designated Orion East by 
Maddalena et al. (1986), this cloud is at a velocity of 
~ lkms -1 and may have little connection to the main 
Orion B cloud, though it seems kinematically related to 
the Orion-Eridanus superbubble. Orion East is discussed 
in the next section. 

Orion B was divided into two sections, and distances 
estimated for each. The area located below b = —14° 
which contains the CO peaks associated with NGC 2023, 
NGC 2024, NGC 2068 and NGC 2071 was designated Re- 
gion 1. The region of more diffuse gas above b = —14° was 
designated Region 2. When the foreground-background 
assignment was made for Region 1 there were three stars 
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Fig. 6. (a) Wco map of Orion B. The integration range is < ulsr < lSkms -1 . Contours start at 1.68 K kms -1 (3a) 
and subsequent levels are at 8, 13, 23, 33, 43, 63, 83, 103, 143 and 183o\ The horizontal dashed line at b — —17° marks 
the nominal boundary between the Orion A and Orion B clouds. The horizontal dotted line at b = —14° separates 
the clouds into two regions for which distances were estimated. Emission from "the Scissors" , marked by the dashed 
box was not included in the velocity-latitude map. (b) Velocity-latitude map of Orion B. The integration range is 
100° < I < 210° and the contours start at 0.514 Karcdeg (3er) with subsequent levels at 8, 13, 23, 33, 43, 53, 63, 73, 
93, 113 and 133er. 



whose extinction was four times below the value that 
would be expected if they were located behind the cloud, 
but ~ 3 times above the extinction predicted by HI alone. 
The dust associated with the molecular gas in the low 
velocity tail of Orion B can probably account for the ex- 
tinction in excess of that expected from HI, suggesting 
that the three slightly reddened stars are located within 
the low velocity tail of Orion B which extends away from 
the main cloud toward the Sun. Two distances to Region 
1 were estimated, a near distance of ~ 420 pc and a far 
distance of ~ 520 pc. The distance to Region 2 was found 
to be ~ 390 pc, though this estimate is rather uncertain as 
only a small number of suitable Hipparcos stars are found 
toward this part of Orion B. 

A kinematic model where the stellar winds from the 
Ori OB lb subgroup have reached Orion B at its low- 
latitude and low-longitude corner is consistent with the 
observed distances. The gas at the interface has been 
compressed laterally and accelerated away from the Sun. 
This gas is located at the far distance. However, a small 
fraction of the gas has been accelerated toward the Sun. 
This gas, which does not appear to be fragmented (c.f the 



2.3 — > 5.5kms _1 channel map - Fig. 0), is estimated to 
have a mass of 1.2 x 10 4 Mq, or about 15% of the total 
mass of the cloud. 

4-- 3. Orion East 

The contours in Fig. [7| are those of velocity- integrated CO 
in the cometary Orion East cloud, located toward the high 
latitude part of Orion B. The Ha emission from the region, 
shown in false colour in Fig. indicates that the ionized 
gas is streaming around the cloud. The pressure of this 
gas may have triggered the recent star-formation that is 
suggested by the presence of at least five TTauri stars 
(Herbig & Rao 1972; Cohen & Kuhi 1979) in the head of 
Orion East. 

Orion East may be at a distance of only 120 pc - signif- 
icantly closer than Orion B and the only previous distance 
estimate of 500 ± 140 pc by Herbst (1982). The new esti- 
mate is based upon the parallaxes of only three Hipparcos 
stars, but since they all had small relative errors and could 
be assigned unambiguously, the new estimate is probably 
robust. In the kinematic model of the Orion-Eridanus bub- 
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Fig. 7. Ionized gas around Orion East (L1622). The 
colour field represents the smoothed continuum subtracted 
Ha emission mapped by the Southern Ha Sky Survey 
(SHASSA, Gaustad et al. 2001). Velocity-integrated CO 
emission from the globules is represented by the superim- 
posed contours. 



ble of Bally et al. (1998), Orion East has been accelerated 
toward the Sun and is on the near side of the bubble. The 
new distance estimate is in agreement with this model. It 
also agrees with estimates concerning the distance to the 
Orion-Eridanus bubble based upon X-ray intensity maps 
(Guo et al. 1995) which indicate that the near side of the 
bubble is at 159 ± 16 pc and the far side is at 226 ± 24 pc 
(Guo et al. 1995) or ~ 400 pc (Reynolds & Ogden 1979). 

Orion East is also noteworthy because it lies beyond 
the projected edge of the Ha emission associated with 
Barnard's Loop and the Orion-Eridanus bubble. Bally 
et al. (1998) argued that the Orion-Eridanus ionization 
boundary had receded, and that Barnard's Loop (and its 
faint Ha extensions into Eridanus) delineate the current 
ionization boundary of the Orion-Eridanus bubble, but in 
the past the influence of the OB association must have 
extended beyond this boundary in order to shape Orion 
East and the L 1617 globules into their present cometary 
form. However, Fig. shows that ionizing gas continues 
to stream around Orion East, and so a more likely expla- 
nation is that Barnard's Loop traces the ionization front 
of the near side of the Orion-Eridanus bubble, which is 
much closer to the Sun than the Orion-Monoceros molec- 
ular clouds. 
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Fig. 8. The Scissors. The colourscale represents the emis- 
sion weighted mean velocity and the contours the velocity- 
integrated CO intensity. The integration range is — 1 < 
^lsr < 16kms _1 and the contours are at 1.8, 3.6, 5.4, 
7.2, 9.0, 10.8 and UTKkms" 1 . The two solid lines mark 
the positions of the position- velocity slices in Fig. The 
cross marks the position of zero offset in these slices. 

4-4- VDB 49: the Scissors 

At I = 202°, b — —14° is an unusual molecular cloud that 
we have dubbed the Scissors. The cloud covers an area of 
^2.2 deg 2 , which at a distance of ~ 150 pc is equivalent to 
~ 15 pc 2 , and it has a mass of ~ 140 M@. The distinctive 
two limb structure and unusual kinematics of the cloud 
are summarized in Fig. [SJ which shows the weighted mean 
velocity field of the cloud as a colour field with contours 
representing the velocity-integrated CO emission superim- 
posed. The solid lines mark the positions and orientations 
of the two position-velocity slices, one along each limb, 
that are presented in Fig. [!|] 

The Scissors appear to be a large cometary cloud with 
a single head that points toward the centre of the OB lb 
subgroup and a forked tail that extends toward the Galac- 
tic plane. However, the two branches of the tail have 
very different kinematics as shown by the position- velocity 
slices, Fig. ED The lower limb (slice 2) has a steep, smooth 
velocity gradient along its entire length. This gradient is 
largely absent in the upper limb (slice 1). The apparent 
velocity discontinuity at the position where the two limbs 
intersect suggests that the Scissors are two separate clouds 
along the same line of sight. However, the probability of 
the chance alignment of two similarly shaped and oriented 
clouds 14° from the Galactic plane is small. 

It is not possible to determine the distance to either 
branch via the constraining star method as there are no 
Hipparcos stars that can be reliably assigned to the fore- 
ground of the Scissors. However, there are two stars, one 
for each branch, that can be assigned to the background 
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Fig. 9. Position-velocity slices through the Scissors. The 
orientation of the slices and the zero offset position 
are shown in Fig. |H1 Each slice is integrated over a 
width of 0?375 wide. The contours in both plots start 
at O.lKarcdeg with subsequent steps at each additional 
O.lKarcdeg. 

and these indicate that neither branch of the Scissors is 
further away than ~ 300 pc. We tentatively adopt a dis- 
tance of ~ 150 pc to both. It is noted that the Scissors does 
not display a strong correlation between Ha emission and 
Wco comparable to that of the nearby cloud Orion-East, 
which has a similar position close to Barnard's Loop and 
a more secure distance. 

A possible explanation for the unusual velocity struc- 
ture of the Scissors is suggested by an examination of 
the kinematics of the surrounding gas. Fig. ^| presents 
a position-velocity map for a wide (~ 6°) strip which in- 
cludes the Scissors, the low longitude tip of the Orion A 
cloud and a number of small cometary clouds. Panel (a) 
shows the velocity-integrated CO map rotated 55° clock- 
wise about the position / = 206°, b = —17°. The inte- 
gration range for the position-velocity plot is shown by 
the two horizontal dashed lines. The positions of OB stars 
from the Ori OB la and lb subgroups (Brown et al 1995) 
are plotted with the same symbols as in Fig.Q] The Wco 
map is presented in this way so that a direct comparison 
with the position- velocity plot, panel (b), can be made. 
The kinematics of these clouds indicate that the Scissors 
forms part of a semi-circular structure that is consistent 
with expansion about the Orion OB la subgroup. The ra- 
dius of the ring is ~ 7°, which corresponds to a radius of 
~ 50 pc at a distance of ~ 400 pc. An expansion velocity 
of 8kms indicates that the ring is ~ 6Myr old, approx- 
imately half the age of the OB la subgroup (Blaauw et al. 
1991). 

In this model of an expanding ring the high velocity 
(~ 12kms _1 ) clouds, characterized by the tip of Orion A, 
represent the most distant part of the ring and the low 
velocity gas (~ — 1 kms -1 ) in the Scissors and the globule 
at x — 8°, y — —5° in Fig. ^J|(a) represents the closest 
part of the ring. This is consistent with the model used 
to explain the kinematics of Orion East, but the actual 
distances derived using the constraining star technique in- 




Fig. 10. Position-velocity slice through Ori OB la. (a) 

Wco map rotated 55° clockwise about the position I = 
206°, b = -17°. The positions of O and B stars from the 
Ori la and Ori lb subgroups are plotted with the same 
symbols as Fig. ^ The two horizontal dashed lines delin- 
eate the position-velocity slice, (b) Position-velocity map 
for the Ori OB la region. The logarithmic contours begin 
at -0.6 (~ 2a) and then proceed in steps of 0.3. 



dicate that the distance between the nearest and furthest 
part of the ring is ~ 300 pc. The distance estimates may 
be in error, as they are based upon only a small number 
of stars, but it is also possible that the ring is older than 
6 Myr and that because of its orientation its actual size is 
underestimated. 

4-5. NGC 2149: a giant expanding ring. 

Between the Mon R2 cloud and the high longitude end of 
the Orion A cloud are a group of molecular clouds associ- 
ated with the reflection nebula NGC 2149 (van den Bergh 
1966). These clouds are at a significantly higher velocity 
(~ 12kms _1 ) than those in Orion A (~ 5 kms -1 ). In most 
projections these clouds appear to be unconnected. How- 
ever, a three dimensional display of the data with velocity 



12 



B. A. Wilson et al.: A Uniform CO Survey of the Molecular Clouds in Orion and Monoceros. 




Fig. 11. (a) Weighted mean velocity field of the NGC 2149 region. The integration range is 1 < wlsr < 16kms _1 and 
the contours are at 1.7 (3cr), 2.8, 7.3, 10.1, 13.0, 15.8, 18.6 and 21.4Kkms~ 1 . The solid ellipse marks the position of 
the expanding ring model fitted to the observations. The solid line shows the orientation of the position-velocity slice 
through the ring that is presented in (b), and the cross marks the location of the projected centre of the ellipse which 
is also the zero offset position. In (b) the slice was integrated over 1?25, the maximum extent of the projected ellipse. 
The contour levels start at 0.2 with subsequent levels at 0.5, 0.9, 1.2, 1.7, 2.3, 2.8, 3.4, 4.0, 4.5, 5.1 and 5.7Karcdeg. 



as a third dimension indicates that the clouds in this re- 
gion may be connected to form a ring. This structure is 
on a much larger scale than the expanding shells located 
within Orion A. 

The structure of the CO emission from NGC 2149 
was modelled as an expanding ring. Fig. ^] presents the 
fit of the model to the observations, (a) is the observed 
peak velocity field of NGC 2149, the ellipse representing 
the position and orientation of the model ring, (b) is a 
position-velocity diagram along the projected major axis 
of the model ring. The two plots in Fig. suggest that 
the model is a reasonable fit to the data. 

The total kinetic energy of the ring is estimated to 
be ^-Mring^exp = 3 x 10 48 ergs, where it is assumed that 
all of the gas in the NGC 2149 region is expanding, i.e. 

M ring = -Mngc2149 = 12 - 3 x 10 3 M©- At a distance of 
400 pc the radius of the ring is w 40 pc, which corresponds 
to an expansion age of ~ 9Myr. The size and energetics 
of the NGC 2149 ring are comparable to those of the A- 
Orionis ring (Lang et al. 2000). However unlike the A- 
Orionis ring, the NGC 2149 ring is not associated with a 
region of excess Ha and there are no stars comparable to 
A-Ori to drive its expansion, which may therefore be the 
result of a supernova. 



4.6. Mon R2 

Mon R2 is a massive star-forming cloud between Orion A 
and the Galactic plane. Its distance is claimed to be 
approximately 800 pc, significantly greater than that of 
Orion A or Orion B, but this estimate is quite uncer- 
tain. The constraining star technique provides only very 
rough upper and lower limits for the distance to the 
cloud, placing it only in the range 600-1100pc. This esti- 
mate is consistent with the photometric determinations of 
830±50 (Racine 1968; Herbst & Racine 1976) and 950 pc 
(Racine & van den Bergh 1970) . Mon R2 subtends an 
area of ~ 15.3 deg 2 (~ 3000 pc 2 ) and has a mass of 
~ 0.75 x 1O 5 M0. Both the size and mass are compara- 
ble to those of Orion A. It is also located at the same 
distance below the Galactic plane, ~ 170 pc, and it has 
a similar velocity gradient that runs roughly parallel to 
the Galactic plane and in the opposite direction to that 
expected from Galactic rotation (Hughes & Baines 1985). 

Mon R2 was first mapped in CO by Kutner & Tucker 
(1975), then by Maddalena et al. (1986) and, at higher 
resolution, by Xie & Goldsmith (1994) with the FCRAO 
14 m telescope. Mon R2 has also been surveyed in 13 CO 
with the AT&T Bell Laboratories 7 m telescope (Miesch 



B. A. Wilson et al.: A Uniform CO Survey of the Molecular Clouds in Orion and Monoceros. 



13 



-9° r • • ' • I ' '-'3" I - " p-» "i ■ ■ ■ ■ 



-10° - 



-12° - 



a -13° 



-14° 




-is- mm 



217° 216° 215° 214° 213" 212° 211° 210 c 
Galactic Longitude 

(a) i pri 1 I.I I I .1 

5 15 25 35 45 




217° 216° 215° 214° 213° 212" 211° 210 c 
Galactic Longitude 

(b) \ f . . ^-ZWH 

8 9 10 11 12 13 



Fig. 12. Mon R2 in CO. (a) Wco map of Mon R2. The integration range is 5 < ulsr < 20kms 1 and the contours 
are at 1.7 (3d), 4.7, 7.6, 10.5, 13.4, 16.3, 22.1, 27.9, 33.7, 39.6 and 51.2Kkms -1 . The two solid lines mark the position 
velocity slices in Fig. US The cross at I = 213?71, b = -12?59 marks the position of the Mon R2 outflow, (b) Mon R2 
emission weighted mean velocity field with Wco contours superimposed. 



& Bally 1994). The high resolution surveys, and a smaller- 
scale multiline study of the high-density central core (Choi 
et al. 2000), have shown that Mon R2 has a unusual in- 
ternal structure where most of the high density gas has 
been cleared from the central core and is constrained to a 
~ 30 pc hemispherical shell. CO observations have also 
revealed the presence of two outflows; a large 200 M Q 
molecular outflow (Loren 1981; Wolf, Lada & Bally 1990; 
Meyers-Rice & Lada 1991) that originates from the main 
core and a much smaller outflow associated with the 
GGD 12-15 core (Little et al. 1990). 

Fig. El(a) shows the integrated Mon R2 CO emission 
from the present survey. Mon R2 apparently contains two 
main structures: most of the emission comes from a pro- 
nounced spur that extends from the central core toward 
higher latitudes and the Galactic plane, or from the low- 
latitude side of the cloud that forms the roughly hemi- 
spherical shell (most clearly seen above the 16 K kms con- 
tour) identified by Choi et al. (2000). Extending along the 
spur are three discrete knots of emission; the largest one, 
which is also the one closest to the main core, contains the 
GGD 12-15 core and associated outflow. 

4.6.1. The Mon R2 shell 

Xie & Goldsmith (1994) reported that a large fraction 
of the molecular gas in Mon R2 is confined to a large 
(~ 30 pc) hemispherical shell, the projected centre of 
which is close to the NGC 2182 reflection nebula. This 



shell is obvious in the two CO position-velocity slices in 
Fig. El because they have been integrated only over the 
position of the core and not the whole cloud. The slices 
go through the position of the main outflow (/ = 213?7, 
b = — 12° 6) and are orientated as shown by the solid lines 
in Fig. 112( a). In Fig. El the main lobe of the outflow ap- 
pears as a distinctive spike. Slice (1) runs roughly per- 
pendicular to the major axis of the outflow, along the 
boundary between the low and high emission sides of 
Mon R2. Fig. El( a ) is dominated by the outflow spike 
but the semi-circular velocity structure that characterises 
the Mon R2 shell is also apparent. Slice (2) runs along the 
major axis of the outflow, as reported by Meyers-Rice & 
Lada (1991); the associated position- velocity plot is shown 
in Fig. 1131 (b). In addition to the outflow spike there is a 
large velocity gradient along the slice which characterises 
the Mon R2 shell. The direction of the velocity gradient 
is not oriented in the same direction as the outflow, which 
accelerates gas at negative offsets toward lower velocities. 
This suggests that the outflow is not responsible for the 
shell. Furthermore, it is clear from Fig. El that the main 
core is a coherent part of the velocity structure of the shell, 
and is therefore likely to be located within it. 

4.6.2. The Mon R2 spur 

Xie & Goldsmith (1994) suggested that the sharp edge of 
the spur was evidence of recent formation, plausibly by 
shock compression. The corrugations that appear along 
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Fig. 13. Two position-velocity slices through Mon R2. 
The positions and orientations of the slices 1 and 2 are 
shown in Fig.ll2l(a). Each slice was integrated over ~ 0?25. 
Contour levels are indicated by the colour bars. 

the spur are also consistent with the kinematic instability 
of shock front development (Elmegreen & Combes 1992). 
Since the main Mon R2 outflow at I = 213?7, b = -12?6 
is the most conspicuous energy source in the cloud, it is 
worth asking if the spur is related to it. 

The spur originates from a position very close to the 
main Mon R2 outflow (marked as a cross in Fig. El 
and it is conceivable that gas entrained by the outflow has 
been aligned along the spur by some mechanism. However, 
for any mechanism to be credible it must also explain two 
significant observations. First, the orientation of the out- 
flow, as determined by Meyers-Rice & Lada (1991), differs 
by ~ 45° from the orientation of the spur, and second, 
the base of the spur is at least 0?5 from the outflow. Jets 
driving smaller outflows have been observed to precess 
(e.g. HH 315 Arce & Goodman 2002), and this mecha- 
nism could explain the observed differences between the 
orientation of the clumps and the outflow, and the offset 
between the base of the spur and the outflow. Neverthe- 
less, the length of the spur (~ 30 pc) may be too large for 
this to be realistic. Wolf, Lada, & Bally (1990) estimated 
the age of the outflow to be 1.5 x 10 5 yr, so for the outflow 
to have swept the most distant clump out to 30 pc the gas 
would have had to have a velocity of 300 km s -1 . Although 
this is comparable to the speed of the jet, it is much faster 
than the speed of the gas entrained by the outflow and so 
it is therefore more likely that another, older, mechanism 
has aligned the spur. 

The projected centre of the shell (I = 213?9, sb = 
— 12? 2) is about halfway between the two dense cores (the 
main core and the GGD 12-15 core) that contain the two 
main Mon R2 outflows. When compared with the centroid 
of the Mon R2 shell the spur appears to be cometary, 
with a well defined rounded head centered on the GGD 12- 
15 core and a tail that extends toward low latitudes. Xie 
& Goldsmith (1994) estimated the dynamical age of the 
shell to be ~ 4 x 10 6 yr and its total kinetic energy to 



be 5-20 x 10 48 ergs. A similar amount of energy deposited 
over about the same length of time might produce the 
observed structure of the spur. 

4.6.3. Star-formation in Mon R2 

Mon R2 is a very active region of star-formation with 
two main condensations centered on the positions of the 
largest bipolar outflows (Shimmins, Clarke & Ekers 1966; 
Beckwith et al 1976; Thronson et al. 1980; Cohen & 
Schwartz 1980; Little, Heaton & Dent 1990). Herbst & 
Racine (1976) suggested that there have been two distinct 
phases of star-formation in Mon R2, the first of which pro- 
duced the many A- and B- type stars and associated reflec- 
tion nebulosity in the cloud. This phase occurred at least 
6 x 10 6 years ago, consistent with the dynamical age of 
the hemispherical shell (Xie & Goldsmith 1994). The sec- 
ond period started 10 4 -10 5 yr ago (Beckwith et al. 1976; 
Loren 1977; Xie & Goldsmith 1994) and it is continuing 
within the dense cores and shock compressed ridges that 
are associated with the bipolar outflows, compact HII re- 
gions, H2O masers and IRAS point sources of the cloud. 
The correlation between the location of star-formation in 
Mon R2 and the dense gas that was shock-compressed by 
the previous generation implies triggered star-formation. 

4-.7. The Southern Filament and the Crossbones 

Extending from close to the Mon R2 cloud to the Galac- 
tic plane is the Southern Filament, a remarkably straight 
and narrow molecular cloud, more than 10° in length 
and generally less than 1° wide. The filament has a rela- 
tively smooth gradient characterized by higher velocities 
at higher latitudes. The velocity gradient can be seen in 
the emission-weighted mean velocity map and position- 
velocity slice, Fig.^l(a) and (b). The location and orien- 
tation of the position-velocity slices is shown by the solid 
line in Fig. OK a )> an d its nominal centre by the cross at 
I = 220°, b= -9°. 

Distances were determined for two parts of the South- 
ern Filament, the regions above and below 6^—7°; both 
found to be at ~ 460 pc. However, these estimates were 
based upon only three and five stars respectively and 
consequently have large relative errors (see Table |2J). De- 
spite the uncertainty in the distance estimate it seems 
more likely that the Southern Filament is related to the 
NGC 2149 ring and high longitude end of Orion A than 
Mon R2. However, there is not a clear boundary between 
the Southern Filament and Mon R2 in either projection 
or position-velocity space. Possible reasons for the kine- 
matic similarity between the gases at different distances 
are explored in §5\ 

The distinctive Crossbones clouds, first identified by 
Maddalena et al. (1986), at I w 219°, b w -10°, may sim- 
ply be the superposition of two clouds along the line of 
sight, but the similar kinematics of the two limbs suggests 



B. A. Wilson et al.: A Uniform CO Survey of the Molecular Clouds in Orion and Monoceros. 



15 




226° 224° 222° 220° 218° 216° 6° 4° 2° 0° -2° 

Galactic Longitude Offset 

(a) flu 



Fig. 14. (a) Emission-weighted mean velocity map of the Southern Filament with Wco contours superimposed. 
Contours start at 1.69 K km s -1 (3cr) with subsequent steps every 1.69Kkms _1 . The solid line shows the orientation 
of the position- velocity slice presented in (b), with the cross marking the position of zero offset, (b) Position- velocity 
plot for the Southern Filament. The data was integrated over a ~ 1° band centered on the filament. Contour levels 
are indicated by the colour bar. 



that this is unlikely. There is some evidence, presented 
in Fig. 1151 that star-formation is currently taking place 
within the Crossbones, near the position where the two 
limbs intersect: (a), the 100 /im infrared emission mea- 
sured by IRAS, shows several knots of emission consis- 
tent with the presence of embedded young stellar objects 
(YSOs). Most promisingly the knot at the intersection of 
the two limbs, I « 219?2, b w — 9?9 is coincident with the 
position of an apparent bipolar outflow suggested by the 
well defined spike in the velocity-latitude map in (c). 

4-8. The Northern Filament 

Fig-ED( a ) is the Wco m &P of a second long filament, com- 
parable to the Southern Filament described above, that is 
located between the high- latitude end of the Orion B cloud 
and the Galactic plane. This cloud was first mapped in CO 
by Maddalena et al. (1986), who called it the Northern Fil- 
ament. However, this name is slightly misleading because 
although the Northern Filament is similar in length and 
orientation to the Southern Filament, it is neither as nar- 
row nor as straight. It is possible that this cloud would not 
have been described as a filament at all had it not been 
aligned parallel to the Southern Filament. 

The distance to the Northern Filament is estimated as 
~ 390 pc by the constraining star technique, the same as 
the distance to the high latitude part of Orion B. In addi- 



tion to being at the same distance, the Northern Filament 
and Orion B have very similar morphologies and kinemat- 
ics along their adjacent edges, suggesting that they were 
once part of a coherent structure. 

The Northern Filament is relatively large and massive: 
it subtends an area of 15.3 deg 2 on the sky and has a mass 
of 0.13 x 10 6 M Q . Despite its size, the cloud is remarkably 
quiescent. Molecular clouds of this size typically have star- 
formation efficiencies of a fraction of one percent, which 
corresponds here to w 200 young stars. There is nothing 
in the Northern Filament to indicate that so much star- 
formation is taking place: there are no T Tauri stars, HII 
regions, reflection nebulae, bipolar outflows or IRAS point 
sources with colours similar to those of EYSOs. 

The kinematics of the Northern Filament, as summa- 
rized by the velocity-latitude diagram, Fig. I16l fb'). are 
fairly unremarkable, which is not surprising considering 
its quiescence. Unlike the Southern Filament, there is no 
large scale velocity gradient, and the whole cloud is con- 
strained to velocities close to ~ lOkrns -1 . Part of the 
reason why the velocity-latitude plot appears so smooth 
is because it has been integrated over a relatively large 
longitude range, 200° < I < 210°; narrower slices do re- 
veal a more complex kinematic structure. In particular 
there are several positions where the small-scale velocity 
structure of the cloud matches that of Orion B at the same 
longitude. 
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Fig. 15. (a) Map of IRAS 100 emission from the Crossbones. (b) Wco map of the Crossbones. The integration 
range is 6.0 < ulsr < 18.0kms _1 . The white (14.5 < v L sr < 18.0kms _1 ) and black (7.0 < u L sr < 9.5kms _1 ) 
contours show the distribution of high and low velocity gas respectively. The two vertical dashed lines show the 
integration range for the velocity-latitude map presented in (c). The crosses mark the positions of colour selected 
IRAS point sources that are potentially EYSOs. The solid line delineates the inferred major axis of the apparent 
bipolar outflow that is centered on / = 219?2, b = — 10?0. (c) Velocity-latitude map of the Crossbones. The integration 
range is 219?2 < I < 219?6kms _1 and the contours start at O.lKarcdeg with subsequent steps every additional 
O.lKarcdeg. The arrows at b = — 10?0 mark the position of the velocity spike associated with the apparent outflow. 



5. Large-scale structure of Orion-Monoceros. 

The similarity between the Orion A and Mon R2 clouds 
has already been noted fSection l4.tjf l. They are at approxi- 
mately the same distance from the plane, they have about 
the same size and mass, and they both have large-scale ve- 
locity gradients in a direction opposite to that of Galactic 
rotation. Fig. ll7l fh) is a position- velocity slice through the 
emission peaks of the Orion A and Mon R2 clouds. The 
orientation of the slice is delineated by the dashed line in 
Fig. El( a ) which shows a composite of two of the Orion- 
Monoceros channel maps (Fig. yjl. The map is presented 
in this way to emphasise that the slice is parallel to the 
Northern and Southern Filaments and that it passes along 
emission spurs in both Orion A and Mon R2 that are also 
parallel to the filaments. The morphological and kinematic 
similarity between Orion A and Mon R2 is striking. Both 
clouds are home to giant molecular outflows, as charac- 
terized by the emission spikes at offsets of -5? 5 and 3°, 
both have lobes of high- velocity gas on their low Galactic 
latitude sides and both have spurs of emission that extend 
from their central core toward the Galactic plane and are 
parallel to the Southern Filament. 



Mon R2 is more that ~ 300 pc further away than 
Orion A and has not usually been considered to be part 
of the Orion complex. However, when confronted by the 
numerous similarities between Orion A and Mon R2 it 
is difficult to avoid speculating that these clouds, and by 
extension all the clouds in the region, share a common 
origin. 

Heiles (1998) proposed that the formation of Orion- 
Monoceros was associated with the expansion of the Vela 
supershell (GSH238+00+09), which may have originated 
from the cluster Collinder 121 (Cr 121) located near 
I ~ 235°, b ~ -10°, u LS R = 17±4kms- 1 (Blaauw 1991), 
and at a distance of 546 ± 30 pc (Hoogerwerf et al. 1997). 
The Vela supershell encloses an enormous irregular bub- 
ble in the region where it has blown out of the Galactic 
plane, one side of which runs approximately parallel to the 
Southern Filament, before becoming more spread out at 
lower latitudes. 

Heiles (1998) claimed that the Vela supershell follows 
the expansion law R oc t 31 , meaning that the walls of the 
supershell had travelled ~ 80% of the way to their present 
location at half the age. Thus the Vela supershell would 
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Fig. 16. (a)Wco map of the Northern Filament. The integration range is 4 < wlsr < 16kms _1 . The contours start 
at 0.97Kkms _1 (3<r) with subsequent steps every additional 3.22Kkms _1 (10<r). (b) Velocity-latitude map of the 
Northern Filament. The integration range is 201° < I < 210° and the contours are at 0.4, 1.2, 1.9, 2.6, 3.3, 4.0, 5.4, 
6.9 and S.SKkms" 1 . 



have overrun the position of the Orion OB association at 
about the time that the OB la subgroup was formed. 

Bally et al. (1998) proposed that a collision of the 
Vela supershell with a fossil remnant of the Lindblad Ring 
may have played a role in triggering the formation of the 
proto Orion-Monoceros complex. They suggest that the 
ISM was first compressed when it was swept up by the 
Lindblad Ring supershell; as the expansion of this super- 
bubble slowed the dense gas in the shell was compressed 
further by the impact of the Vela supershell, enhancing 
the tendency of the gas to collapse via gravitational insta- 
bility and form stars. 

The structure of Orion-Monoceros close to the 
Orion OB association is now dominated by the local effects 
of stellar winds, supernova explosions and the expansion 
of the Orion-Eridanus bubble, but much of the underly- 
ing large-scale structure is consistent with the formation 
sequence described above. In particular the alignment of 
the Southern Filament and the Orion A and Mon R2 spurs 
with the wall of the Vela superbubble does not appear to 
be purely coincidental. 

6. Summary and Conclusions 

This paper presents the results of a new CO survey of the 
Orion-Monoceros complex carried out with the Harvard- 
Smithsonian 1.2 m telescope. The new observations are 



more sensitive and have better sampling than the only 
previous large survey of the region (Maddalena et al. 
1986). Distances to the majority of the clouds in Orion- 
Monoceros are determined using the constraining star 
technique and the Hipparcos catalogue. The new distances 
aided the three dimensional analysis of the molecular gas 
and it relationship to the other phases of the ISM. The 
main conclusions of this analysis are as follows: 

1. Stellar winds and ionizing radiation from the massive 
stars of the Ori OB lb and OB Id subgroups are re- 
sponsible for changes to the morphology and kinemat- 
ics of the molecular gas in Orion A and Orion B today. 

2. The kinematics and three dimensional structure of 
Orion A and Orion B are consistent with stellar wind 
driven compression centered on Ori OB lb. Knots of 
compressed gas are coincident with the most active 
star-forming regions in Orion A and Orion B, suggest- 
ing that triggering may have occurred. 

3. The Orion-Eridanus superbubble of hot Ha emitting 
gas inflated by the Orion OB association and de- 
lineated in projection by Barnard's loop has over- 
run many small clumps of dense gas. These globules 
have cometary shapes caused by their interaction with 
the expanding bubble, and many display signs of en- 
hanced (triggered) star-formation in their heads. How- 
ever, only one globule, Orion-East, displays a fairly 
unambiguous connection to the ionizing radiation of 
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Fig. 17. (a) Composite Wco map of Orion-Monoceros. 
The area within the box marked by the dotted line has 
been integrated over the range 12.0 < ulsr < 15.3 km s -1 , 
with the remainder of the map integrated over the range 
8.8 < i>lsr < 12.0kms _1 . The contours are the same as 
the channel maps in Fig. The dashed line delineates 
the position-velocity slice presented in (b) and the cross 
marks the position of zero offset, (b) Position- velocity 
slice through Orion A and Mon R2. The contours start 
at -0.9 (~ 3a) and go up in steps of 0.3 (log Karcdeg). 



Barnard's Loop, and this globule is significantly closer 
than the main molecular clouds of the complex. There- 
fore, we conclude that Barnard's Loop delineates the 
closest part of the Orion-Eridanus bubble and not the 
part adjacent to Orion-Monoceros. 
4. The most distant cloud in the complex, Mon R2, is 
quite similar to the large molecular cloud associated 
with the Orion nebula that is significantly closer. Both 
clouds have spurs of emission that are aligned with the 
Southern Filament (and the Vela supershell). These 



features suggest that the clouds of the region share a 
common origin. 
5. The new CO survey supports the proposals by Bally et 
al. (1998) and Heiles (1998) that the formation of the 
Orion-Monoceros complex was triggered by the pas- 
sage of the Vela supershell and possibly its interaction 
with a preexisting supershell associated with the Lind- 
blad Ring. 
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